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Abstract 

The generally accepted methods of predicting heat- 
transfer coefficients often give values which are in poor 
agreement with the measured coefficients for single 
phase hydrogen flowing through straight symmetrically 
heated circular tubes for the extreme conditions en- 
countered in the NERVA or Phoebus-2 nozzle. A 
thorough study of the 4600 experimental heat-transfer 
coefficients for hydrogen flowing through straight tubes 
is made which results in an equation which predicts co- 
efficients with acceptable accuracy over a wide range of 
conditions. The prediction equation for straight tubes is 
modified to predict coefficients for hydrogen on both the 
" swept concave surface" and the "upswept convex sur- 
face" of curved symmetrically heated circular tubes. 
Without further modification the equation predicts heat- 
transfer coefficients which are in good agreement with 
values measured in asymmetrically heated noncircular 
channels simulating the throat region of the NERVA and 
Phoebus-2 nozzles. These tests were run under the ac- 
tual anticipated heat flux, pressure, temperature, and 
flow rate. The prediction equations are recommended 
for use in predicting heat -transfer coefficients in the 
cooling passages of the actual NERVA and Phoebus-2 
rocket nozzles. 


Introduction 

The extreme conditions encountered in regenera- 
tively cooled nuclear rocket nozzles produce severe 
heat-transfer problems in the coolant passages. An ef- 
fective method of predicting heat-transfer coefficients in 
the cooling passages is essential to the optimization of 
any nozzle design. Of particular concern is the high 
heat flux throat region where fluxes of 20 Btu/(sec)(in. 2 ) 
and higher are reached. 

A number of investigations have been conducted with 
single phase hydrogen flowing turbulently through tubes. 
The range of conditions in these Investigations have ap- 
proximated those encountered in the cooling passages of 
a nuclear rocket nozzle. These studies resulted in sev- 
eral correlation equations each of which was limited to a 
particular range of conditions. A recent study found 
that a single equation can be used to predict heat- 
transfer coefficients over a much greater range of con- 
ditions than was previously possible. f 1 ) 

In this presentation the most recent straight tube 
prediction equation(l) is first compared with the two 
prediction equations most widely used in rocket nozzle 
design. 12, 13) Th e straight tube prediction equa- 
tion! 1 ) is then modified to include the effects of curva- 
ture and compared to some existing experimental data 


for single curved tubes. Finally the suggested applica- 
tion of these prediction equations to the cooling passages 
of a nuclear rocket nozzle is discussed. 

Heat-Transfer Coefficients in Straight Tubes 

A study of all available turbulent single phase hydro- 
gen heat-transfer data has been reported in Ref. 1. 

These 4622 data points resulted from 10 investigations 
using symmetrically heated straight circular tubes with 
a straight unheated approach length. ( 2 to u ) The com- 
plete range of experimental conditions covered by the 
investigations is shown in table I. The study reported in 
Ref. 1 gave as a result the correlation equation 


Nu b = 0. 023 Reg- 8 Prg- 4 ^jexp - (o. 57 - (1) 

which was tested for the range of inlet pressure and in- 
let temperature shown in Fig. 1 . The use of the 
temperature-entropy diagram is a convenient method of 
showing the location of measured inlet pressure and tem- 
perature in relation to the saturation lines and critical 
pressure p c and critical temperature T c . 

The hydrogen data are separated into regions 1 and 
2 as shown in Fig. 1. These regions are the result of 
the study of the effect of inlet temperature Tj, inlet 
pressure p i( and the transposed critical temperature 
(the temperature at which the specific heat at constant 
pressure reaches a maximum) T*. In region 1, 87 per- 
cent of the 3674 heat -transfer coefficients predicted by 
equation (1) deviates less than ±25 percent from the 
measured values. Region 2 is defined by 45° R < Tj <T* 

p c < pj < 530 psia and Is often referred to as the 
near-critical region. In region 2, 40 percent of the 948 
predicted heat -transfer coefficients deviates less than 
±25 percent from the measured valves. At present there 
is considerable doubt about the transport properties in 
this region. 


There are two equations commonly used in pre- 
dicting heat -transfer coefficients in the cooling passages 
of a regeneratively cooled nuclear rocket nozzle. One 
of these is the Hess and Kunz equation^ 2 ) 


NUf 


: 0.0208 Re° -8 Pr 


f' 4 (l + 0-01457 ^ 


( 2 ) 


The other one Is the modified Hess and Kunz eous- 
tion( 10 and 13) H 


Nu f = 0. 0208 C L Re^' 8 Pr l 


■f' 4 ( 1 + 0 -01457 ^j 


(3) 
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The quantity C L varies nonlinear ly as follows: 


Coolant 

C T 

temperature. 


°R 


50 

2.0 

55 

1.73 

60 

! 1.48 

65 

1.26 

70 

1.07 

75 

.93 

80 

.87 

85 

.85 




7 

0.05 


( 4 ) 


for the concave side. 


The convex side of a curved tube is rarely given any 
consideration, but in this study the reciprocal of the Ito 
factor was used to modify equation (1) to give 


N % 


0.023 Re£’ 8 Pr; 


O.4/M / 

b \T h J exp ~ 57 " 

* h(»)] 


1.59 

x/D 

-0.05 


L ) 


( 5 ) 


for the convex side. 


It°0 r lant temperatures above 85 ° R C L is constant 

Equations (1) to (3) are used to predict heat- 
transfer coefficients which can be compared to the most 
recent experimental data for the range of conditions en- 
countered in nuclear rocket nozzles. The ratio of these 
predicted coefficients to the measured values are shown 
as a function of the ratio of wall to fluid bulk tempera- 
tures in Fig. 2. Figure 2 shows that 94 percent of the 
heat-transfer coefficients predicted using equation (1) 
deviated less than ±20 percent from the measured values 
compared to 42 percent using equation (2) and 66 percent 
using equation (3). ^ 

Heat-Transfer Coeffi cients in Curved Tubes 

t ,, ^ eat transfer measurements have been reported for 
both the concave or "swept" side and the convex or 
unswept" side of curved circular tubes. ( 14 > Equa 
tion (1) was first tried without any change from the form 
that correlated straight tube data to predict heat- 
transfer coefficients for both the concave and convex 
sides of the symmetrically heated test sections The 
ratio of the heat-transfer coefficient predicted by equa- 
tion (1) to the measured coefficient is shown as a func- 
tion of temperature ratio in Fig. 3. As might be expec- 
ted, the predicted heat-transfer coefficient is lower than 
e measured value on the concave side and higher than 
e measured value on the convex side. The calculated 
heat-transfer coefficient is as small as one-half the ex- 
perimenta 1 value on the concave side and as large as al- 
most twice the experimental value on the convex side. 
These effects of curvature on the heat-transfer coeffi- 
c tents have been noted elsewhere 

„„„ Eq . Uatl ° n (1) ls modi fied in this paper with the Ito 
correction factor for curvature* 17 ) to give the equation 

N« b - 0.023 


. The resul ts of using equations (4) and (5) to predict 
the heat-transfer coefficients on the concave and convex 
sides of a tube are shown in Figs. 4(a) and (b), respec- 
tively. The experimental data In Fig. 4 is the same data 
used in Fig. 3. Approximately 90 percent of the pre- 
dicted heat-transfer coefficients deviated less than 
±20 percent from the experimental values. 

It appears that the Ito correction does improve the 
correlation and that equations (4) and (5) will predict 
heat-transfer coefficients on the concave and convex 
sides of curved tubes, respectively. 

Thus far, both the straight and curved tube data 
were for symmetrically heated circular tubes. In a 
rocket nozzle the cooling passages are noncircular and 
they heated from only one side. Some data are available 
or asymmetrically heated circular and noncircular pas- 
sages curved to simulate the throat and high flux region 
of the NERVA and Phoebus-2 nozzles. < 10 > Heat-trlfer 
coefficients predicted by equation (4) were compared to 
these experimental values for conditions as near as pos- 
sible to those for the actual nozzle. 

Figure 5 shows the ratio of the heat-transfer coeffi- 
cients calculated by equation (4) to the experimental co- 
efficients as a function of the dimensionless distance 
from the entrances, m Fig. 5(a) the ratio of coefficients 
are for ft e NERVA contour and in Fig. 5(b) the ratio is 
for the Phoebus-2 contour. The agreement between pre- 
dicted and measured heat-transfer coefficients is con- 

^ ^ b jL better . f ° r thC Rhoebus-2 configuration than for 
the NERVA configuration. The ratio of coefficients 
seems to vary more erratically for NERVA data than for 
the Phoebus -2 data. 

One possible explanation for this behavior might be 
that for the NERVA contour. The calculated temperature 
drop through the wall is up to ten times the temperature 
difference between the inside wall and the hydrogen For 
this condition an error of 2 percent in the measured out- 
side wall temperature can result in a 30 percent error in 
the heat-transfer coefficient. A 2 percent error in out- 
side waU temperature will result in only 10 percent in 
the Phoebus-2 contour where wall temperature drop is 
never more than three times the temperature difference 
between the inside wall and the hydrogen. Another ex- 
planation might be due to the difference in curvature of 
the two contours. 

An important use of the heat-transfer coefficient is 
m the prediction of wall temperatures. Equation (4) is 
use to predict the inside wall temperatures shown in 
g. 6. The experimental wall temperatures and the 
wall temperatures predicted by a modified Hess and 
Kunz equation are taken directly from Ref. 10. Only the 


2 



two runs shown in Figs. 6 (a) and (b) were shown in 
Ref. 10 and, therefore, are the only ones shown in this 
paper. As shown both equations usually predict wall 
temperatures higher than the experimental values with 
the modified Hess and Kunz values being more conser- 
vative. 


feet of curvature would diminish with x/D rather than 
change abruptly at the point of tangency. The curvature 
effect appears to be present at the last instrumented sta- 
tion shown In Fig. 6 (b) which is 2 diameters downstream 
of the tangency point. 


Heat-Transf er Coefficients in Rocket Nozzle 
Cooling Passages 

Regeneratively cooled rocket nozzles are made up 
of noncircular passages formed to give the desired area 
ratios, and are essentially combinations of straight and 
curved tubes. The success demonstrated in predicting 
heat transfer coefficients in straight with tubes using 
equation ( 1 ) and curved tubes using equations ( 4 ) and ( 5 ) 
encourages the use of these equations in predicting heat 
transfer coefficients in the coolant passages of a rocket 
nozzle. 


To predict the heat transfer coefficients in the cool- 
ing passage of a nozzle, the following equations are rec- 
ommended: 


Equations ( 1 ). (4), (5), and ( 6 ) are the best available 
for the prediction of heat-transfer coefficients in the 
cooling passages of hydrogen cooled rocket nozzles. The 
use of equations (1), (4), (5). and ( 6 ) should allow the 
design of more efficiently and reliably cooled rocket 
nozzles. 


Symbols 

c p specific heat of gas at constant pressure 
D inside diameter of test section 


mass flow rate per unit cross-sectional area 
local heat -transfer coefficient 


Entrance of Coolant Passag e 
Straight entrance x/D > 1 

""t, ■ °- K ’< 4 (|f)™P - («.B - $£) ,1) 

45° and 9 0° angle bend entrance x/D > 5 
N % .0.023KeJ- 8 p I »-^^.( 0 . 5 ,.l 5 |) 

X ('* F 1?) I®) 

where Fj = 3.5 for the 45° angle bend and 5 for the 90° 
an gIe bend entrance. (18 1 Fl differs from the value 
given' > by 1 . 4 (the value for F x for a straight en- 
trance and which is included in the exponent of T w /T|j). 


k thermal conductivity of gas 

Nu Nusselt number, hDA 

Pr Prandtl number, c p p/k 

P absolute static pressure 

q rate of heat transfer to gas per unit area 

R radius of curvature 

Re b bulk Reynolds number, GD/p b 

Re f modified film Reynolds number, PfV b D/p f 

r inside radius of passage 


Throat Section (Concave Curvature) 



Exit Section (Convex Curvature) 

Nu,, = 0. 023 Re°- X- 4 (^)exp - (0. 57 - 


X 



Any Straight Sections 

NU b = °‘ 023 Re b 8pr b 4 " ( 0> 57 ' Vd) 


(4) 

(5) 
( 1 ) 


At present it is not known how far the effects of cur- 
vature on the heat -transfer coefficient will extend down- 
stream of the point of tangency where a curved tube be- 
comes straight. It is reasonable to assume that the ef- 


T temperature 

T transposed critical temperature (temperature at 

which specific heat of fluid at constant pressure 
reaches a maximum) 

v velocity 

x distance from entrance of test section 

p absolute viscosity of gas 

v kinematic viscosity of gas, p/p 

P density of gas 

Subscripts: 

b bulk (when applied to properties, indicates evalu- 

ation at bulk temperature T b ) 

c critical 


3 


f film (when applied to properties Indicates evaluation 
at film temperature Tf) 

i inlet 

w wall 
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x/D 

2.0 

252 

T w/T b 

1.1 

27.6 

T\, o R 

4.5 

573 

Pj, psia 

18 

2500 

q, Btu/(sec)(in. 2 ) 

0.036 

27.6 

Re b 

7500 

13 800 000 

T w , °R 

53 

5600 


TABLE I. - RANGE OF EXPERIMEN- 
TAL CONDITIONS 
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Figure 1. - Range of hydrogen inlet temperature and inlet pres- 
sure for which equation (1) has been experimentally checked. 
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(a) h Calculated using equation (1) (Ref. 1). 
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(b) h Calculated using equation (2) (Ref. 12). 
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V T b 

(c) h Calculated using equation (3) (Refs 10 
and 13). 

Figure 2. - Variation of the ratio of calculated 
to measured heat transfer coefficients with 
temperature ratio. Straight tube data from 
reference 10, x/D from 6. 7 to 33. 9, heat 
flux from 6. 4 to 27. 6 Btu/(sec)(in. h 






E-3436 




(a) Concave surface, equation (1) 
(ref. 1). 



V T b 


(b) Convex surface, equation (1) (ref. 1). 


Figure 3. - Variation of the ratio of calculated 
to measured heat transfer coefficients with 
temperature ratio. Curved tube data from 
reference 14, symmetric heating. 
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(b) Convex surface, equation (5). Present study. 

Figure 4. - Variation of the ratio of calculated 
to measured heat transfer coefficients with 
temperature ratio. Curved tube data from 
reference 14, symmetric heating. 
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(b) Phoebus-2 contour, h calculated using equation (4) (present 
study. ' 

Figure 5. - Variation of the ratio of calculated to measured heat 
transfer coefficients with axial distance from entrance. 

Curved tube data from reference 10, asymmetric heating. 



(a) NERVA contour. 



0 25 30 35 40 45 50 55 60 65 70 


x/D 

(b) Phoebus-2 contour. 

Figure 6. - Variation of inside wall temperature with axial dis- 
tance from entrance. Data from reference 10. 


NASA- CLEVELAND, OHIO E-3436 



